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The analytic form of the lattice quark propagator is used to derive the functional form for short distance mesonic 
correlators. These are then used to calculate "Continuum Model" Ansatze which comprise of a pole, representing 
the ground state, plus a contribution for the excited states, coming from the short distance behaviour. These are 
compared to Monte Carlo data. 



1. Introduction 

On the lattice, the usual method for studying 
hadronic physics is to determine the two-point 
correlation functions of hadronic operators. 
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where M is the mass of the ground state, and 
r is channel dependent. The presence of the de- 
caying exponential has led to investigations being 
focused on the large Euclidean time "tails" of the 
two-point functions, and data at small times be- 
ing discarded. 

The question was asked, "Is there a way of ob- 
taining information from short times?". The an- 
swer was "Yes" - a QCD Sum Rules approach. 
Lattice QCD allows the determination of correla- 
tion functions deep within the non-perturbative 
region. QCD Sum Rules (QCDSR) can be utiHsed 
in the near-perturbative regime where the excited 
states cannot be ignored. Both are techniques 
which have been widely employed to deepen the 
understanding of strongly interacting systems, 
and have been combined in the work of Leinwe- 
ber[0j|], and Allton and Capitani 

In this study, we extend the work of ||l|-^ by 
applying the analysis with a continuum limit ex- 
pansion of the lattice quark propagator, derived 
ini. 

The lattice data used is UKQCD generated, 
and is (quenched) Wilson at (3 — 6.0 with a vol- 
ume of 16^ X 48. 



2. QCDSR Continuum Model 

In this section we introduce the QCDSR-CM 
and apply it to the continuum quark propagator. 
The QCDSR-CM replaces a discrete spectrum of 
excited states with a continuous spectrum of spec- 
tral density p, and introduces a sharp threshold 
in the energy scale, sq, marking the onset of the 
continuum. 

The basic object studied is the continuum OPE 
for QCD quark propagator: 
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where m is the quark mass. This is substi- 
tuted into the (Wick-contacted) 2-point function 
at zero-momentum: 

G2{tf''''=Jd'xTr{S,{xfP''rS,{-xfP''r}.i3) 

It is useful to express the time-sliced correlation 
function, G2{t)'^^^, in the spectral representa- 
tion: 
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and the OPE spectral density is calculated by 
means of an inverse Laplace transform of G2^^ . 
We invoke the QCDSR-CM by setting the thresh- 
old So in the energy scale, so that the excited 
states' contributions to G2{t) is given by only the 
energies above that scale, 

p{s) = Z S{s - M) +ie{s~ so) p(s)°^^. (5) 

We have to include the ground state in order to 
obtain the full correlation function. This is in- 
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eluded as a (5 function in p{s) in Eq.(|^). So: 



G2{t) = Ze-*" + e / Pisr'^'' e'^' ds. (6) 

There are four parameters in the fitting ansatz: 
Z — > normaUsation of ground state 
M — > ground state mass 

^ — > normalisation of excited states, introduced 
to allow for lattice distortions ||^ 
So — > continuum threshold. 
The OPE expansions for the degenerate mesons 
and baryons using the continuum quark propaga- 
tor are listed in The analysis was extended 
to non-degenerate mesons in |5j. 

3. Wilson Quark Propagator 

In terms of loop-divergent behaviours, a coor- 
dinate space formulation of a lattice theory, has 
the benefit over its momentum analogue of high- 
lighting the short-distance phenomena which are 
the origin of the poles. In the continuum, the an- 
alytic expression for a fermion propagator in po- 
sition space is well known, Eq.(|2|). On the lattice 
however, the standard representation involves in- 
tegrals over Bessel functions, and proved difficult 
to analyse in the continuum limit. Liischer and 
Weisz successfully analysed the massless propa- 
gator in 4-dimensions While in the mass 
dependence of the propagator at x = was de- 
rived. However, in Q|, the asymptotic expansion 
for the modified Bessel function of the first kind, 
/„, which appears in the expression for the lattice 
propagator, (using Schwinger's integral represen- 
tation). 
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was calculated, leading to the following expres- 
sion for the Wilson lattice propagator: 



(a;2 
' 2ram*Ko 



2 



m(x 



2ram* — — - 



2\2 



m(x ) 



-I- 4m^(l + ram)'^^^^K2 



2^ 2 



m[x ) 



(8) 



where Ki and K2 are modified Bessel functions of 
the second kind. This expression was then used in 
a QCDSR-CM style approach following the pro- 
cedures outlined in section 2. 
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Figure 1. Effective mass plot for the 2- 
exponential and 0{a) QCDSR-CM fits for the 
vector channel a,t k — 0.1530. 
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Fit Z M xVd.o.f 

Type 

Pseudoscalar 

K = 0.1550 

1- exp 1.81(6) X 10-2 0.2986(15) - 

2- exp 1.91(4)xl0-2 0.3014(12) 10/15 
0{a) 1.7(3)xl0-2 0.294(3) 8/15 

K ^ 0.1540 

1- exp 2.12(5) X 10-2 0.3645(11) - 

2- exp 2.28(4) X 10-2 0.3685(10) 24/15 
0(a) 1.6(14) X 10-2 0.356(14) 18/15 

K = 0.1530 

1- exp 2.46(5) X 10-2 0.4236(9) 

2- exp 2.67(3) X 10-2 0.4281(8) 42/15 
0(a) 2.1(4)xl0-2 0.416(3) 30/15 

Vector 

K = 0.1550 

1- exp 6.7(5) X 10-3 0.426(4) 

2- exp 8.5(2) X 10-3 0.440(2) 160/15 
0(a) 5.9(2) X 10-3 0.419(3) 30/15 

K = 0.1540 

1- exp 8.6(3) X 10-3 0.467(3) 

2- exp 1.03(2) X 10-2 0.479(2) 210/15 
0(a) 7.5(2) X 10-3 0.461(2) 37/15 

K = 0.1530 

1- exp 1.05(3) X 10-2 0.509(2) 

2- exp 1.25(2) X 10-2 0.5193(13) 240/15 
0(a) 9.3(2)xl0-3 0.5020(15) 41/15 



Table 1 

Values of Z, M and for the fitting ansatze. 



4. Wilson Correlators 

In this section, the work of 0-|3) and 0| is com- 
bined by applying a QCDSR-CM approach to the 
continuum limit expansion of the Wilson lattice 
propagator, i.e. we use S^^^ = in Eq.(^). 

To obtain a closed expression for G2 we first ex- 
pand the modified Bessel functions for short dis- 
tances and then perform the spatial integral. This 
procedure was performed for all mesonic chan- 
nels, with general (non-degenerate) constituent 
quarks, and for the momenta, p = and p =- 
(1, 1, 1). We obtain, e.g., for the zero-momentum 
pseudoscalar channel the following expression: 

Gl^itT^^ = -^(l-t-2ram) 
3^2 

Values for the Z and M fitting parameters to- 
gether with the for the pseudoscalar and vec- 
tor channels are given in Table |ll The QCDSR- 
CM fits for the Wilson correlators are labelled 
by 0(a). The fitting window for the 1-exp fit 
was t — [12, 20] and for the 2-exp and 0(a) fits 
was t — [2,20]. Note that the 0(a) fits perform 
better than the 2-exp fit in every case, and that 
the improvement is more significant for the vector 
channel particularly as the mass increases. 

5. Conclusions 

In this work we have derived lattice mesonic 
correlation functions for the Wilson action, in- 
cluding 0(a) effects. Our Ansatze fit lattice 
Monte Carlo data better than conventional "2- 
exponential" fits. Future work includes fitting of 
heavy mesons, and deriving analogous formulae 
for the Clover action. 
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